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Abstract 

Recently, the collisionless expansion of spherical nanoplasmas has been analyzed with a new 
ergodic model, clarifying the transition from hydro dynamicdike to Coulomb-explosion regimes, 
and providing accurate laws for the relevant features of the phenomenon. A complete derivation of 
the model is here presented. The important issue of the self-consistent initial conditions is addressed 
by analyzing the initial charging transient due to the electron expansion, in the approximation of 
immobile ions. A comparison among different kinetic models for the expansion is presented, showing 
that the ergodic model provides a simplified description, which retains the essential information on 
the electron distribution, in particular, the energy spectrum. Results are presented for a wide range 
of initial conditions (determined from a single dimensionless parameter), in excellent agreement 
with calculations from the exact Vlasov-Poisson theory, thus providing a complete and detailed 
characterization of all the stages of the expansion. 
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I. INTRODUCTION 



The irradiation of solid targets with ultraintense lasers can induce the prompt formation 
of hot dense plasmas, which rapidly expand into a vacuum, as predicted by J. Dawson in 
1964 1]. In the case of planar targets (employed for ion acceleration [2j), the physics of 
the expansion has been extensively studied under a variety of conditions, using different 
analytical and numerical approaches 3), 0, 0, O, tJ. In contrast, the expansion of spherical 
plasmas [such as the nm- to /im-sized plasmas generated upon interaction of ultraintense 
lasers with atomic or molecular clusters (cf. Q, (9), Q, 11, 12])] have not been analyzed as 
thoroughly. A deep knowledge of the expansion (accounting for the self-consistent dynamics 
of ions and electrons) can be relevant in particular situations where accurate control over the 
expansion is necessary, examples being the double-pump irradiation of deuterium clusters 
aimed at tailoring the ion dynamics so as to induce intracluster fusion reactions [13( , or the 
bio molecular imaging with ultrashort X-ray pulses jjjj], where expansion control is needed 
to avoid significant damages of the sample before the typical imaging time. In fact, accurate 
solutions for spherical expansions exist only for ideal cases, such as the Coulomb explosion 
(CE) [jjj] of a pure ion plasma, which occurs when all the electrons are suddenly swept 
away from the cluster by the laser field. In opposite conditions, when most of the electrons 
are heated by the laser but not stripped from the cluster, hydrodynamic models can be 
employed 16j to estimate the basic features of the expansion; in the quasineutral limit, a 
kinetic solution for the adiabatic expansion of plasma bunches into a vacuum has also been 
derived 17|. However, in more general situations (for example, in experiments with large 
clusters, containing millions of atoms), a significant violation of charge neutrality occurs, 
even though a relevant fraction of the electrons remains bound to the cluster. In such 
conditions (corresponding to plasma radii on the order of the electron Debye length, or 
less), the expansion process is strongly dependent on the self-consistent dynamics of ions 
and trapped electrons and it can be described accurately only by kinetic models, based on 
the Vlasov-Poisson (VP) theory. A remarkable example of numerical solution of the VP 
equations for the expansion problem in a spherical geometry can be found in Ref. ?]|, for 
the particular case in which the motion of both ions and electrons is purely radial. 

Recently, the Authors presented a kinetic analysis of the collisionless expansion of spher- 
ical plasmas driven by hot electrons, based on a peculiar ergodic model, which accounts 
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for the radial motion of the ions and for the three-dimensional motion of nonrelativistic 
electrons [3]. In the present paper, the model is derived in detail, and its validity is tested 
against reference solutions of the full VP equations (here obtained using ad-hoc numerical 
techniques). Furthermore, a procedure to determine the self-consistent initial conditions for 
the expansion within the framework of the model is presented. As an accurate knowledge 
of the initial space-charge distribution is fundamental to describe correcly the long-term 
plasma expansion, the initial charging transient, during which the faster electrons leave the 
cluster core, is analyzed resorting to different models, thus providing deeper physical insights 
and validating the technique. 

The results presented here provide a complete characterization of the expansion dynamics 
of spherical nanoplasmas, which can be useful in the interpretation of recent experiments 
with clusters, either irradiated with intense IR lasers or with VUV/X-ray sources pjj], where 



conditions may be far from those of a pure CE (lOj. In particular, the different behavior 
observed in the ion energy spectrum for different values of the electron temperature allows a 
clear identification of the transition from hydrodynamic-like to CE expansion regimes, thus 
defining the range of validity of the CE approximation. 

II. KINETIC MODELS FOR THE EXPANSION 

In the electrostatic, nonrelativistic limit, the dynamics of a collisionless ion-electron 
plasma is described rigorously by the Vlasov-Poisson (VP) set of equations: 
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where / e (r, v, i) and /i(r, v, t) are the distribution functions in phase space for electrons 
(having mass m and charge — e) and ions (having mass M and charge Ze), respectively, 
and $ (r, t) is the electrostatic potential (<& is set to zero at infinity, so that the energy of a 
single electron, e = |mt> 2 — e$, is negative if it is trapped). In the following, the attention is 
focused on the expansion of a plasma sphere (with initial radius Rq) composed of cold ions 
(with initial uniform density n i0 ) and hot electrons [with initial uniform density n cQ = Zn i0 
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and arbitrary energy distribution p e {e)]. The general initial conditions for Eqs. ([I]) can be 
cast in the form 

fco (r, v) = n e0 g (v ) 9 ( 1 - ] 

/io(r,v) = ™io<M v ) ( 1 " ^ ) 

where G is the Heaviside step function and g is an arbitrary function of v, such that p e {e) 
= 167rVoi?o/3 (2e/m 3 ) 1/2 #[(2e/m) 1/2 ]. In the present paper, for simplicity, only the ideal 
situation (commonly adopted in the literature 3, 0, S, Q, E]) of an initially neutral plasma 
with Maxwellian electrons is considered, in which all the information on the electron heating 
by the laser pulse is contained in the initial electron temperature, T . The function g in Eq. 
(J2J) is then defined as g (v) = [m/ (27r/cBTo)] 3 ^ 2 exp [— mv 2 / (2/cbT )]. As can be readily proved 
by writing Eqs. ([I]) and (|2J) in nondimensional form, in this case the dynamics of the system is 
fully determined by the dimensionless parameters Zm/M and T = Z^B^o/ecE = 3A D0 /.Rq, 
being Ado the initial Debye length for the electrons, and ecE = ^f e2 R^ n to the maximum ion 
energy attainable in the case of pure CE for a sphere of ions. With the use of the initial 
conditions ([2]) the electrons are supposed to be instantaneously heated by an infinitely short 
laser pulse, without expanding. However, in principle, any initial space-energy density for the 
electrons could be employed (for example, linear superpositions of Maxwellian distributions, 
as well as initially nonneutral distributions), as resulting from a realistic model of laser- 
matter interaction. Nonetheless, the use of the reference initial conditions ([2]) provides a 
simple way to obtain useful physical insights on the influence of the combined effects of 
variations of electron energy and cluster features. 

The expansion process is split in two stages: first, a rapid expansion of the electrons, 
which leads to a VP equilibrium before the ions move appreciably; second, a slow expansion 
of the plasma bulk, driven by the positive charge buildup formed in the first stage. Due to 
the large mass disparity between ions and electrons, a simplified model can be derived, in 
which the ions are assumed as immobile during the former process, whereas the electrons 
can be considered as instantaneously at equilibrium with the electrostatic potential during 
the latter stage. A self-consistent theoretical framework can be developed, which allows one 
to determine accurately both the initial equilibrium and the bulk expansion, by treating 
the electron dynamics as a sequence of equilibrium configurations (with frozen ions). The 
model is obtained by exploiting the functional relation existing between n e and $ at equilib- 
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rium, and by calculating the energy variation of the electrons under the hypothesis of slow 
variations of $ in time. 

Equilibrium solutions of the Vlasov equation for the electrons must depend on r and v 
only through the invariants of motion. Since a spherical symmetry has been assumed and the 
electrostatic force is central, the only invariants of motion to be considered are the Hamil- 
tonian 3tf (r, v) = mv 2 /2 — e$ (r) and the angular momentum, L = mr x v. Consequently, 
the equilibrium distribution function can be written as f e (r, v) = F (3f (r, v) , mr x v). If 
a generic space point r = re r and a generic velocity v = v r e r + v j_e_|_ are considered, the 
phase-space density is given by F y— (v 2 + v\) — e$ (r) , Jzf , e r x e_i_J being Jz? = mrv± the 
absolute value of L. Due to the spherical symmetry of the system (and, in particular, the 
symmetry with respect to any rotation with respect to e r ), the phase-space density cannot 
depend upon e_i_, and, consequently, F depends only on 3f and Jzf. 

The energy-angular momentum distribution, a c (e,£), can be defined as 

a c (e,£)= Jj F [34? (r, v) , ^ (r, v)] 5 [3? (r, v) - e] 5 [jgf (r, v) - i\ drdv = 
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where R\ (e, I) and R 2 (e, £) (R\ < R2) are the radial turning points, i.e., the values of r such 



that e — 
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e<3> (r) = 0. The quantity a c (e,£) AeAi represents the number of electrons 



having energy in [e, e + Ae] and absolute value of the angular momentum in [£, I + A£]. The 
electron density, n e , can be written as 



n e (r) = ^JJa c (e,£)V (r, e, £; {$}) ded£, 
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is such that V(r;e,£)Ar gives the probability, for an electron with energy e and angular 
momentum £, to be found in [r, r + Ar]. If time variations of due to the ion motion, are 
slow with respect to the period of the of the radial oscillation of the electrons, the mean 
value of de/dt can be evaluated as 
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i.e., by using the ensemble average of -— . This is equivalent to preserve the value of the 
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Equations ©-(EJ), coupled with Poisson's equation and Newton's equation for the radial 
motion of the cold ions, provide a self-consistent model for the collisionless expansion of a 
finite-size plasma in the case of spherical symmetry. 

In kinetic theory, there is a precise relationship between time scales and the proper 
number of parameters to be used to describe correctly a given phenomenon: in the case 
of the plasma expansion, the VP system (CQ) for f e and /; allows one to follow precisely 
the expansion dynamics on the time scale of the fastest particles; to study the ion ex- 
pansion, a quasi-equilibrium model, Eqs. ([3])- ([6]), can be used, in which the stationary 
solution of the Vlasov equation for the electrons is employed. In fact, as the Vlasov model 
is noncollisional, it does not contain a physical mechanism leading towards the equilib- 
rium (the equations are time-reversible). To justify the use of the equilibrium distribution 
/ (r, v) = / [Jf (r, v) , Jzf (r, v)], one must suppose that the stationary solution of Vlasov 
equation is a good representation of the real electron distribution, once high-frequency fluc- 
tuations are eliminated; formally, this can be performed by introducing a dissipation mech- 
anism (i.e., a suitable collision term into the Vlasov equation). In general, an approximate 
kinetic model can be regarded as the result of introducing a particular collision term. For 
example, by using a binary collision term with sufficiently high collision frequency, f e tends 
towards the Maxwell-Boltzmann distribution [i.e., / (r, v) = Const • exp (— Jt? (r, v) /k&T)], 
in which all the information is restricted to the temperature; in this case, a hydrodynamic 
description is obtained, whose domain of validity is confined to situations where T < 1. 
For larger values of T , the use of a proper energy spectrum is fundamental; in fact, the 
energy distribution presents a cutoff for e = (for e > the electrons are not confined and 
their stationary density must vanish) and this fact is hardly compatible with a Maxwellian 
distribution having a non-negligible fraction of electrons with e > 0. 



Within this framework, the approach of Ref. 18] can be introduced by considering a 
collision term of the form 

J (f e ) = -v (/ e - / e ) , / e = i- I U (r, vh, t) df2 (8) 



where f2 is a unit vector and v represents the collision frequency; the specific value of v is 
irrelevant, as long as \ jv is much smaller than the characteristic time of the ion expansion. In 
this case, the kinetic Vlasov equation for the electrons is replaced by the collisional equation 

dt dr m dr dv 

-vf + ^ff(r,vn,tjdn, (9) 

where v is the collision frequency and Q is a unit vector. The collisions do not alter the 
electron energy, but change randomly their direction, driving f e towards an equilibrium 
distribution having the form f e (r, v) = / [J>tf (r, v)], a sort of ergodic density such that each 
electron has an equal probability to be found in every point of the hypersurface of phase- 
space having equation J4? (r, v) = e. This is different from the usual ergodic distribution 
of the statistical mechanics, in which the state of a system of N particles can be found 
with equal probability on the hypersurface J$? (r 1; r 2 , r N ; v 1; v 2 , v N ) = Const of the 
complete, 6iV-dimensional phase-space. In fact, the time derivative of the electron entropy 
S c = — jj f e log (/ e ) drdv can be written in the form 

dt 

which is always nonnegative unless / e = / e . Therefore, a necessary condition for the distribu- 
tion to be stationary is that f e must not depend on f2 . Finally, the equilibrium distribution 
is a function of and 5£ that does not depend on fl, and, consequently, it is a function of 
Jf? only. 

In the following, the approach will be referred as single-particle ergodic (SPE) method. 
According to this approach, the equilibrium distribution function can be written simply as 
f e (r, v) = / [J$? (r, v)], the dependence on Jzf being lost, and Eqs. ©-([5]) are replaced by 

P°(z)=jJ f ( r > v )] 6 ( r ' v ) " £ ] drdv = 
16tt 2 V2 
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where 31 (e) is the integration domain, such that r G @ (e) =>- e + e$ (r) > {for monotonic 
potentials, ^ (e) = [0, R (e)], where e + e<& (i? (e)) = 0}, and the adiabatic invariant, Eq ([7j), 
is replaced by the ergodic invariant [20I 

J" [e (t) , t] = Const ■ / [e + e$(r,t)]^r 2 dr, (14) 

J 9(e) 

defined as the volume of the region of M 6 enclosed by the hypersurface of equation ^mv 2 — 
e$(r) = e. 

As shown in Ref. 18|, the SPE approach provides excellent results for the expansion of a 
spherical plasma in a wide range of the parameter T . There is a number of reasons to explain 
its success, even though it is not easy to quantify their relative importance. First, even 
though hydrodynamic models can provide a qualitative agreement with the real expansion 
dynamics, the SPE model is extremely more flexible in describing the energy distribution of 
the electrons. In addition, in the cases considered here, the initial phase-space distribution, 
Eq. ([2]), is assumed to be an SPE function. Moreover, it must be noticed that the angular 
momentum is invariant only in the case of perfect spherical symmetry. In practical situations, 
perturbations to that symmetry (e.g., an initial shape which is not perfectly spherical, or 
collisions with heavy particles) would cause a mixing in Jzf distribution, and their effect 
could be taken into account by introducing a collision term such as the one in Eq. ([8]). 



III. SINGLE-PARTICLE ERGODIC MODEL 



Under the hypothesis of SPE distribution, a self-consistent model for the expansion of 
a spherical plasma can be formulated, starting from Eqs. (|TT]) -( ri4|) . as follows [18]. A 
Lagrangian approach can be used both for the ions (which move in radial direction, starting 
from the initial position r , with zero velocity) and for the electrons (whose energy e evolves 
in time starting from the initial value eo), by determining the ions trajectories n (ro,t), the 
electron energies e (eo,£), the ion density n ; (r, t), the electron density n e (r,t), the electron 
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energy distributions p e (e, t), and the potential <3> (r, t) according to the set of equations 
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de/de 
J(e {e ,t),t)=Q, 



where the evolution equations for the radial coordinates of the ions and the electron energies 
are coupled via the Poisson's equation. The expansion dynamics is determined once the 
initial ion density n i0 and the electron energy distribution p e0 are given. In Eqs. (fT5l) . 
the electron density n c is expressed as the sum of the number of electrons having energy 
in [e, e + de] [i.e., p e (e) de = p Cj0 (e )de )], multiplied by the probability for an electron with 
energy e to be found at the radius r, according to the ergodic distribution. For simplicity, 
the ion density n; is written under the hypothesis of no ion overtaking (dri/dro ^ 0}_ 21]; 



however, the model can be easily generalized to include many-branched shock shells 



13 



21 



34 



221 ] and different ion species. 



The set of equations (FIBTj describes the expansion dynamics on the ion time scale; there- 
fore, its numerical solution is much faster than solutions of the full VP model (where the 
electron time scale must be followed). The model is solved by calculating the radial tra- 
jectories of a set of representative ions and the energy variations of a set of computational 
particles. Each computational particle represents a given number of electrons, whose radial 
distribution is given by Eq. (TT3"]) . This description of the energy dependence corresponds to 
a suitable discretization of the integral in Eq. f[T2"|) (which is similar to the description of 
the spatial dependence commonly adopted in particle- in-cell approach 3]). 



IV. CHARGING TRANSIENT 



Since the initial configuration considered here [cf., Eqs. ([2])] is far from equilibrium, the 
proper p Ei o, to be used in Eqs. f|T5l) . must be determined as the equilibrium configuration 
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following the initial charging transient. Apparently, the SPE method cannot be of help 
to this purpose, since it is valid only for sufficiently smooth variations of $ in time, a 
condition which is not met in the early stage, when the hot electrons are suddenly allowed 
to expand (as if a rigid wall, initially confining them, were instantaneously brought to 
infinity). However, a procedure has been envisaged, which makes these equations suitable 
also for the analysis of the initial electron equilibrium, thus allowing the study of the whole 
expansion process (initial charging transient and bulk expansion) within the same theoretical 
framework. Before describing this procedure, the initial equilibrium is analyzed in detail, 
accounting for the full electron dynamics (VP model). 



A. Reference solutions of the collisionless and collisional models 



Due to the importance of an accurate knowledge of the initial equilibrium configuration 
of the electrons for a correct analysis of the plasma expansion, reference results for the 
transient leading to the initial space-charge distribution of the plasma has been determined 
by solving Eqs. (pQ) numerically, in the hypothesis of immobile ions. In the numerical scheme 
adopted, computational particles representative of a given number of electrons are moved in 
space, under the action of the sum of the self-consistent electric field and the electric field 
due to the ion distribution. By resorting to the spherical symmetry of the system, the field 
generated by the electrons is evaluated using Gauss' law, as if each particle were actually a 
spherical shell, thus avoiding the use of a computational grid for solving Poisson's equation, 
and allowing for an infinite radial domain (similar techniques have been used to investigate 
the VP dynamics of one-dimensional plasmas 241 ]: here, the validity of the method has been 
checked through comparisons with reference results from the 3D particle-in-cell code OSIRIS 
25l|). The same framework has been used also to investigate numerically the effect of the 
presence of the collision term, Eq. ([9]), which forces the system towards a SPE distribution. 
Such perturbations are introduced in the model by scattering randomly the computational 
particles, without changing their energy, according to the collision frequency v. 

Figure [1] shows the evolution of the electronic charge contained within the ion sphere, for 
the representative low-temperature (To = 7.2 x 10~ 3 ) and high-temperature (To = 7.2 x 1CT 2 ) 
cases of Ref. [3], as obtained with the pure VP model [y = 0) and with the collisional 
model, Eq. ([9]), (using v = u> pe ). In the collisionless case, the charge transient exhibits small- 
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amplitude oscillations (the simulation parameters have been carefully checked to ensure that 
the oscillations are not due to numerical noise). In the collisional model, for v > u> pe , the 
oscillations are strongly damped and the system rapidly reaches an equilibrium configuration, 
as predicted theoretically. 

B. Method of the barrier 

In order to build a self-consistent ergodic model for the whole expansion process (thus 
avoiding the use of different models to deal with the initial stage), a procedure has been 
devised to determine the equilibrium distribution that follows the initial electron expansion, 
using the same theoretical framework of Eqs. (|TT]) - (fl4|) . To this purpose, the charging 
transient described by the full VP model is replaced by a virtual charging transient, in 
which an external potential barrier, initially confining the electrons, is gradually moved 
from i?b = Rq to infinity with a series of small radial displacements. Each time the barrier 
is moved farther by 5Rb, the new self-consistent potential $ is calculated and the energy of 
the electrons is updated. In order to actually simulate an expansion into vacuum (which the 
real transient is), the electrons and the expanding barrier must not exchange energy, i.e., 
the electron energy must vary only because of $ variations. This implies that the ergodic 
invariant ( fl4l) is not conserved during the initial stage. In fact, should one conserve J when 
displacing the barrier from a given radius i?b to Rb + SR^, the corresponding electron energy 
variation, 5e, would be 

rRb 

5e = -e 5$Q(r,e;{$})dr-5W, (16) 
Jo 

where 6W, defined as 

8W = H [ e + e $ (R h )} Q (R h ; e) 5R h . (17) 

represents the expansion work, done by an electron having energy e, against the expanding 
barrier. Thus, conserving J would cause the overestimation of the electron cooling as the 
system would lose an extra amount of energy corresponding to the expansion work. In 
order to obtain an energy balance equivalent to that of a vacuum expansion, the energy loss 
associated to the expansion work is set to zero in Eq. 0161) . 

The physical process simulated with the barrier method can be thought as an infinitely 
slow expansion during which some external energy source exactly compensate for the ex- 
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pansion work SW against the barrier, or, alternatively, as a series of instantaneous, small, 
displacements of the barrier, where, after each displacement, one waits for a new equilibrium 
configuration to establish. 



C. Drift-diffusion approximation 



As an alternative to solving Eq. (Q, one can consider the drift-diffusion equation 

2 1 d 
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dt dt de Zmv r 2 dr 
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dr 2 dr 
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obtained from Eq. (Q by approximating f e as fo(r,v) + v • fi(r, v) |26[. In Eq. (fT8l) , the 
quantity ^ (r, e, £) represents the space-energy distribution of the electrons (i.e., \I/ (r, e, £) Ae 
is the particle density for electrons with energy in the range [e, e + Ae]). The self-consistent 
potential <3? is determined by solving Poisson's equation 
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Asymptotically, for t — ► oo, the solution approaches a stationary solution of Eq. ffl8|) . vl/^, 
such that 
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dr ar 
By solving Eq. ([2"0]) with respect to Vl/oo, one finds 
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p QO (e) [e-fegoo(r)]^ 
4tt / r /2 [e + e$ 00 (r / )]3dr' 

J 9(e) 



(20) 



(21) 



(poo(e) = 4n f \E' 00 (r, e)r 2 dr is the energy distribution), which corresponds to the SPE dis- 
tribution expressed by Eqs. (JTT]) - (JT3]) . 



D. Results and comparison between models 

Examples of initial equilibrium are now presented and discussed, first referring to the 
two cases of Fig. [TJ then examining the full T -dependence of the principal equilibrium 
parameters. A comparison of the self-consistent equilibrium configuration of the electrons 
after the initial charging transient is made between the exact VP model, Eq. (pQ) , the barrier 
method, and the drift-diffusion model [Eqs. f|T8l) - f|T9l ]. 
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In Fig. [2, the electron density is plotted, along with the corresponding electric field: 
the positive charge buildup at the ion front, AQ, is (a) 12.5% and (b) 38% of the total 
ionic charge eN . Figure [3] shows the equilibrium energy distribution p e , to be used as 
initial condition for the bulk expansion. The excellent agreement between different models 
confirms the validity of the barrier method. Figure H] shows the asymptotic solution of Eqs. 
( TLB"]) and (|T9l) . ^/^(r, e): the corresponding electron density and energy distribution, plotted 
in Figs. [2] and [3], have been calculated as J ^/ OQ (r,e)de and Am f ^/^(r, e)r 2 dr, respectively. 

The dependence of the initial equilibrium on T has been analyzed using both the barrier 
method and the drift- diffusion approximation, for To varying in the range [10 _3 ,1]: the 
equilibrium values of AQ and of the mean kinetic energy of the trapped electrons, S, are 
displayed in Figs. [5] and respectively, along with the corresponding fit laws (obtained in 
Ref. jl8f] using the SPE model), 

£ 

where F^x) = x/(l + s 1 ') 1 ^, and where the coefficient \/6/e in Eq. (122 1) provides a match 
with the analytical results for the planar case 0, ^] in the limit To <C 1. Again, an excellent 
agreement is found between different calculations. 



^.60 (yeWj , (22) 

1- J- 3 .35 (l.86f 1/2 V (23) 



V. BULK EXPANSION 



The self-consistent expansion of ions and electrons has been investigated for a wide range 
of the parameter To by solving Eqs. (fT5j) . having used the barrier method to determine 
the initial equilibrium distribution of electrons. The results of the study reveal that the 
expansion dynamics changes smoothly from a hydrodynamic-like regime (in which the outer 
ions expand first and a rarefaction front propagates inward) to a CE-like regime (in which all 
ions start expanding at the same time), when going from T <C 1 to T ~ 1. Nonetheless, a 
qualitative change in the ion energy spectrum is observed for To — 0.5, marking the transition 
towards a CE behavior. Following the organization of Sec. IIVDI the bulk expansion is first 
analyzed in detail for the two reference cases [cases (a) and (b) henceforth], in which (a) 
T = 7.2 x 10~ 3 and (b) T = 7.2 x 10~ 2 , and then the dependence of the most relevant 
expansion features on T is examined. 
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The evolution of the ion phase-space profile and of the electron and ion densities (starting 
from the initial equilibrium of Fig. [2]), are shown in Figs. [TJandlH], respectively. In case (a) 
(Figs. [T^l and[Sk), ^ ne i° n expansion starts from the periphery and a rarefaction front is 
clearly observed to propagate inward until it reaches the center of the distribution; during 
the expansion, the plasma remains approximately neutral, apart from the thin double-layer 
at the ion front. These features, typical of quasineutral, hydrodynamic expansions, are lost 
in case (b), in which all the ions are promptly involved in the expansion (Fig. [7b) and the 
distribution remains nonneutral during the whole process (Fig. [Hb). In both the 
ions expand, and gain kinetic energy, the electrons cool down and the charge buildup within 
the ion sphere decreases, as illustrated in Fig. [case (a)] and Fig. [TU] [case (b)]. Asymptot- 
ically, the sphere enveloped by the expanding ion front encloses all trapped electrons, and a 
ballistic regime is reached for both species [7j. The self-consistent behavior of the electrons 
strongly affects the ion dynamics and their resulting energy spectrum. In fact, starting from 
the equation of motion of the ions [the first of Eqs. ( fT5|) ]. the asymptotic energy eoo of an 
ion can be written as 

6oc(r ) g(fb,0) | r 1 d q ( ri (r ,t),t) d ^ ^ 



Ze r J r;(r ,t) dt 

where q(r,t) is the net charge buildup enclosed by a sphere of radius r at time t. The first 
term on the right-hand side of Eq. (I24"l) is the ion potential energy, whereas the integral 
term (vanishing for a CE) accounts for the energy loss due to the decreasing charge buildup 
experienced by the ions along their trajectory. Figure [TT] illustrates the evolution of the 
ion energy spectrum towards its asymptotic form: in both cases, the spectrum develops a 
well-defined local maximum far from the cutoff energy. Since this feature is absent in CEs, 
where the asymptotic spectrum is always monotonic (it behaves as e 1 / 2 up to the cutoff 
energy €ce), the maximum in the spectrum is expected to disappear when increasing T 
further. This transition from nonmonotonic to monotonic ion spectra occurs about T = 0.5 
(cf. Fig. fT2l) and marks the transition towards a CE-like behavior. In this sense, To = 0.5 
can be considered as a lower bound for the validity of the CE model. The dependence of the 
maximum (cutoff) ion energy e max on T is shown in Fig. [13j along with the energy value of 
the local maximum in the spectrum, e pea k. The behavior of e max is accurately described by 
the fit law 

Cmax = ^1.43 ( 2.28 Tq 3 ^ 4 J ecE, (25) 
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[T belongs to the same class of functions used in Eqs. (1221) and (123|) ). whereas e pea k exhibits 
the power-law behavior e p0 ak = 0.3TQ- 9 ecE, for T < 0.5. These fit laws can be used to 
provide useful estimates of the initial electron temperature and, hence, of the expansion 
regime. This can be important for the interpretation of experimental ion-spectrum data. In 
fact, for expansion conditions far from a CE, the nonmonotonic behavior of the single-cluster 
ion spectra could affect the total (i.e., arising from all expanding clusters) energy spectrum 
measured in experiments: inparticular, it could lead to nonmonotonic energy spectra, such 
as those presented in Ref. 10J, also for narrow distributions of cluster radii. 



VI. CONCLUSIONS 



The results presented in the paper prove that the collisionless expansion of spherical 
plasmas driven by hot electrons can be analyzed accurately with a kinetic model that de- 
scribes the electron distribution as a sequence of ergodic equilibrium configurations. The 
self-consistent equilibrium that is established after the initial, sudden expansion of the elec- 
trons has been investigated in detail, in the frozen-ion approximation. This equilibrium can 
be determined with great accuracy by replacing the real (fast) transient with an appropriate 
virtual (slow) process, finding excellent agreement with reference solutions of the full VP 
model. This guarantees a highly-precise description of the whole process, thus providing an 
effective tool for the analysis of the expansion dynamics. In particular, a transition in the 
behavior of the ion energy spectrum, when approaching the Coulomb-explosion regime, has 
been identified, and accurate fit laws for the general properties of the expansion, which are 
valid for any value of dimensionless electron temperature (provided that relativistic effects 
are negligible), have been determined. These laws can furnish useful estimates for the inter- 
pretation of experimental data, in particular concerning possible influences of single-cluster 
effects on measured ion spectra. 

Finally, the ergodic model presented here can be readily employed to study more general 
physical situations, such as expansions driven by initially non-Maxwellian electrons. Fur- 
thermore, the model could be extended so as to include relativistic velocities and to account 
for the effects of non- instantaneous electron heating by ultraintense laser pulses. 
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VII. FIGURES 
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FIG. 1: Evolution of the electronic charge contained within the ion sphere (r < Rq), for (a) 
To = 7.2 x 10~ 3 and (b) To = 7.2 x 10~ 2 . Thick gray lines refer to the collisionless case, thick black 
lines refer to the collisional case (y ~ u; pe ). Thin horizontal lines indicate the results obtained 
using the barrier method described in Sec. IIVB1 Units are normalized to the quantities indicated 
in parentheses. 
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FIG. 2: Equilibrium electron density (black, left axis) and radial electric field (gray, right axis), 
for (a) To = 7.2 x 10~ 3 and (b) To = 7.2 x 10~ 2 . Solid lines refer to results from the ergodic model, 
markers to results from the full VP model, and dotted lines to results from the drift-diffusion model 
of Eqs. (|18p and (|19p (in the plots, the curves obtained with the SPE model and those obtained 
with the drift-diffusion model are undistinguishable without magnification). Units are normalized 
to the quantities indicated in parentheses. 
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FIG. 3: Equilibrium energy spectrum of trapped (e = mv 2 /2 — e<3? < 0) electrons, for (a) Tq = 
7.2 x 10~ 3 and (b) To = 7.2 x 10~ 2 . Solid lines refer to results from the ergodic model, markers to 
results from the full VP model, and dotted lines to results from the drift-diffusion model. Units 
are normalized to the quantities indicated in parentheses. 
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FIG. 4: (Color online) Equilibrium distribution in the r — e phase space, ^, as obtained with the 
drift-diffusion model of Eqs. (HI]) and CE|), for (a) T = 7.2 x 10~ 3 and (b) f = 7.2 x 1(T 2 . 
The spatial-energetic distribution is normalized to its maximum value $ mffi and isolevel curves are 
plotted at lO"™/ 2 * 

max (where n is an integer). Units are normalized to the quantities indicated in 

parentheses. 
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FIG. 5: Equilibrium charge buildup as a function of To. Markers are used for results obtained with 
the SPE model (circles) and the drift-diffusion model (asterisks); the solid line represent the fit 
law of Eq. (|22p . whereas the dotted line shows the corresponding power-law behavior for To <C 1. 
Units are normalized to the quantities indicated in parentheses. 
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FIG. 6: Equilibrium total kinetic energyas a function of To. Markers are used for results obtained 
with the SPE model (circles) and the drift-diffusion model (asterisks); the solid line represent the 
fit law of Eq. (|23p . whereas the dotted line shows the corresponding power-law behavior for To <C 1. 
Units are normalized to the quantities indicated in parentheses. 
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FIG. 7: Evolution of the ion phase-space profile, for (a) f = 7.2 x 1(T 3 and (b) % = 7.2 x 1(T 2 . 
Lines refer to results from the ergodic model, markers to results from the full VP model. Insets 
show the evolution of the ion phase-space profile during the early stage of the expansion. Units 
are normalized to the quantities indicated in parentheses. 
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FIG. 8: Evolution of ion (solid line) and electron (dashed curve) density, for (case a, left panels) 
Tq = 7.2 x 10~ 3 and (case b, right panels) T = 7.2 x 10~ 2 . Units are normalized to the quantities 
indicated in parentheses. 
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FIG. 9: Radial trajectory of the ion- front (top panel), evolution of the number of electrons enclosed 
by the ion front (bottom panel, gray), and evolution of the mean kinetic energy of trapped electrons 
(bottom panel, black), for Tq = 7.2 x 1CT 3 . Units are normalized to the quantities indicated in 
parentheses. 




FIG. 10: Radial trajectory of the ion- front (top panel), evolution of the number of electrons enclosed 
by the ion front (bottom panel, gray), and evolution of the mean kinetic energy of trapped electrons 
(bottom panel, black), for To = 7.2 x 1CU 2 . Units are normalized to the quantities indicated in 
parentheses. 
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FIG. 11: Evolution of the ion energy spectrum (from light gray to black), for (a) To = 7.2 x 10 
and (b) T = 7.2 x 10 -2 . Units are normalized to the quantities indicated in parentheses. 




FIG. 12: Asymptotic ion energy spectra for different values of To, compared with the theoretical 
asymptotic spectrum for the CE case (gray curve). Units are normalized to the quantities indicated 
in parentheses. 
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FIG. 13: Cutoff ion energy (black) and location of the maximum in the ion energy spectrum (gray) 
as functions of Tq: circles refer to the SPE model, solid lines to the fit laws in the text. The dashed 
line represents the power-law behavior of e max for To <C 1. Units are normalized to the quantities 
indicated in parentheses. 
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